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Abstract 

This paper is devoted to investigate the cold plasma wave proper- 



f^ ' ties. The analysis has been restricted to the neighborhood of the pair 

C^^ ■ production region of the Kerr magnetosphere. The Fourier analyzed 

^L , general relativistic niagnetohydrodynamical equations are dealt under 

. • I special circumstances and dispersion relations are obtained. We find 

f^ ' the x-component of the complex wave vector numerically. The corre- 

^D . sponding components of the propagation vector, attenuation vector, 

phase and group velocities are shown in graphs. The direction and 

dispersion of waves are investigated. 
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1 Introduction 

Einstein's general theory of relativity describes gravity as a curvature of 
spacetime caused by the presence of matter. Most massive physical objects 
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like black holes distort the geometry of spacetime by their immense gravity. 
The accretion disk around a rotating black hole (most physical in nature) is 
always filled with charged plasma due to the pair production in this region. 
This charged plasma creates an external magnetic field around the black hole 
(the region containing this field is termed as magnetosphere in literature) and 
is governed by the general relativistic magnetohydrodynamical (GRMHD) 
equations. The intrinsic angular momentum of the black hole intends the 
surrounding plasma to act as a magneto-rotational fiuid. The black hole 
spin modifies the rotation of the plasma which fiows in the form of an MHD 
wind. When this wind is perturbed analytically, GRMHD waves are formed. 
These waves propagate to transmit information in the magnetosphere. They 
also help us to detect whether the energy extraction is allowed by the plasma 
state in the magnetosphere. 

There is a large body of literature available [l]-[9] about the perturbations 
and GRMHD waves in the regime of the Schwarzschild black hole. The 
study of rotating black holes has been an interesting subject since they lie 
in the hearts of active galaxies [TU]. Blandford and Znajek |TT] discussed 
the possibility of an electromagnetically driven wind from the rotating black 
hole. Takahashi et al. [12] discussed the MHD infiow and the extraction of 
energy in the form of jets. 

In general relativity, different observers generally experience distinct time 
measures. In 1962, an attempt was made by Arnowitt et al. [TB] who quan- 
tized the gravitational field by a specific foliation (now called ADM 3-1-1 
formalism). Thorne and Macdonald ^^-^ISj used the ADM 3+1 split to 
bridge the gap between the black hole electrodynamics and that of fiat space 
electrodynamics. The formalism was taken over by Holcomb and Tajima 
[TB] and Dettmann et al. [T7] to discuss the quantum description of plasma 
including electromagnetic waves in FRW metric. Buzzi et al. [TS] used 
this access to investigate the wave phenomenon near the Schwazschild event 
horizon by using two component plasma. Recently, Sharif and Sheikh |19j - 
[22] worked on the properties of cold and isothermal plasmas living in the 
neighborhood of the Schwarzschild event horizon. Komissarov [23] discussed 
Blandford-Znajek monopole solution for black hole electrodynamics using 
3+1 formalism. Khanna [24J derived the generalized Ohm's law in 3+1 split. 

The energy can be extracted from the ergosphere via a process first pro- 
posed by Penrose [25]: a particle entering the ergosphere can split in two in 
such a way that one fragment falls into the hole, but the other leaves the ergo- 
sphere with more energy than the original particle. The extra energy comes 



from the hole itself. Leiter and Kafatos [25] discussed the Penrose pair pro- 
duction in massive Kerr black hole's ergosphere. Kafatos [27] discussed the 
gamma ray observations from Penrose powered black holes. Koide et al. [28] 
modeled the GRMHD behavior of plasma flowing into the rotating black hole 
in a magnetic field and showed (using numerical simulations) that energy of 
the spinning black hole can be extracted magnetically. Another way of ex- 
tracting black hole's spin energy is the form of super-radiance. Kokkotas and 
Schmidt [29] discussed quasi normal modes for black holes and relativistic 
stars. They have also discussed the w-waves or gravitational wave modes in 
these scenarios. Furuhashi and Nambu [30] discussed the super-radiance in 
charged Kerr spacetime by numerically solving the Klein- Gordon equation 
as an eigen value problem. Zeldovich [31] considered the waves scattering 
by a rotating conducting cylinder. They proved that Kerr geometry displays 
super-radiance. Zhang [32]- [33] formulated the black hole theory for station- 
ary symmetric GRMHD with its applications in Kerr geometry. He showed 
(taking the specific values of the angular frequency and the x-component of 
the wave vector) that the cold plasma allows the outflow of energy flux. Re- 
cently, Sharif and Sheikh [31]-[SZ] extended this work to cold and isothermal 
plasmas by investigating the real wave numbers. They verified the results 
given by Zhang [33]. This work partially helps to answer the question of 
energy extraction. 

In this paper, we shall discuss the cold plasma wave properties by cal- 
culating complex wave numbers. This provides an alternative to verify the 
results obtained in [37]. We have considered the Fourier analyzed perturbed 
GRMHD equations for the planar analogue of the Kerr metric. These equa- 
tions lead to the dispersion relation and consequently the x-component of the 
complex wave vector for the cold plasma. This complex x-component can be 
split into x-components of the propagation factor (real part) and attenuation 
factor (imaginary part). Further, we evaluate the x-components of the phase 
and group velocities etc. near the pair production region. These quantities 
would help us to understand the modes of dispersion. 

The layout of this paper is as follows. Section 2 provides the description 
of spacetime near the pair production region with some assumptions. Section 
3 consists of the background flow assumptions and related quantities used in 
specification and interpretation of results. In Section 4, we Fourier analyze 
the GRMHD equations for the cold plasma. Section 5 is devoted to the 
numerical solutions of the dispersion relation. It also includes figures and 
related discussion. The observable quantities are given which help us to judge 



the waves properties more efficiently. We shall conclude this discussion in 
Section 6. 



2 Description of Model Spacetime 

The most general line element in 3+1 formalism can be written as 

ds^ = -a^dt^ + -iij{dx' + (3'dt){dx^ + (3^dt), (2.1) 

where a is the lapse function, /3* are the components of shift vector and 7jj 
are the components of spatial metric. All these quantities are functions of 
time and space coordinates. 

We consider the planar analogue of the Kerr metric [31] 

ds'^ = -dt^ + {dx + P{z)dtf + dy"^ + dz^ . (2.2) 

The directions z, x and y are the analogue of the Kerr's radial r, axial (^ and 
poloidal 6 directions respectively whereas t represents the time coordinate. 
Utilizing coordinate freedom, the lapse function is set to unity without loss 
of generality. This has been done to avoid complications which arise due to 
the presence of horizon and hence, the red-shifts. The rotation of the black 
hole imposes a shift vector which is obvious from Eq. (l2.2p . The value of 
the shift function /3 (analogue to the Kerr-type gravitomagnetic potential) 
decreases monotonically from (2; — )■ 00) to some constant value {z — )• —00). 
Equation 02. 2p shows that our assumed /3 is I3{z) = /3(2;)ex which derives 
an MHD wind and extract translational energy analogous to the rotational 
energy of the Kerr metric. The plasma particles, created at 2; = (pair 
production region), are then driven up to relativistic velocities by magnetic- 
gravitomagnetic coupling, as they flow off to "infinity" [z = +00) and down 
towards the "event horizon" {z = —00). 

Our scenario is based on the pair production region and its surrounding 
plasma. Thus we are interested in the magnetospheric plasmas rotating in 
the 2-D fashion i.e. in (zx-plane). The shift is restricted to be dependent on 
only one dimension. It may depend on the other dimensions which obviously 
make the analysis tough but more general. The shift is the effect of black 
hole rotation, of course. The plasma is being produced in the pair production 
region (electron-positron pair production holds and at z = in our scenario) 
outflows from the region towards the event horizon which lies at 2; = —00 and 



towards the outer end of the magnetosphere (the end away from the event 
horizon). The plasma produced in the pair production region falls out from 
the region. This outflow disturbs the neighboring magnetospheric plasma 
in the form of waves. The region z = is not a physical singularity but 
the region where particle-antiparticle pairs are produced from the vacuum 
energy. Since the out flowing plasma is moving in the x2;-plane, so it will 
absolutely disturb the outer plasma in the same plane. 

This spacetime requires a fiducial reference observer (FIDO) (analogous 
to the zero angular momentum observer of the 3+1 split of the Kerr space- 
time), the one with four- velocity perpendicular to the hypersurfaces of con- 
stant time t. The geometrized units will be used throughout the paper. 

3 Plasma Flow and Relative Assumptions 

We consider a rotating background which is assumed to be filled with cold 
plasma admitting the following equation 



fj, = — = constant, (3.1) 

Po 

where /i, po and p are the specific enthalpy, rest and moving mass densities 
of the fluid respectively. This shows that the plasma has vanishing thermal 
pressure and thermal energy. Thus the cold plasma particles are moving 
along timelike geodesies of the spacetime (12.21) . 

The rotating fluid's four-velocity and magnetic field measured by our 
FIDO can be described by spatial vector fields lying in x^-plane 

V = V{z)e^ + u{z)e,, B = B{\{z)e^ + e,}, (3.2) 

where S is a constant. 

For the stationary symmetric background, the relationship between x and 
z components of the velocity vector for the spacetime (12.21) can be written 
as [33] 

V = C + \u, (3.3) 

where C = /3 + Vp is the normalized shift and Vp is an integration constant. 
Further, the mass conservation law in three-dimensional hypersurface is 

apo'-yu = A (constant) , (3.4) 



where 7 = ,^ \ , ,^ is the Lorentz factor. 

' vi— « — V 

We take the assumptions Vp = 1, (3 = tanh(z) — 1 for simphfications. It 
is important to note that our assumed /3 is continuous across z = (i.e., the 
pair production region). This shift may not be vahd for the near horizon. The 
expansion of shift vector /3 imphes that this vector is directly proportional to 
the dimensionless angular momentum J = -^ of the black hole and inversely 
proportional to the radius of the black hole. Thus the spin of the black hole 
can be expressed by the following relation p^j 

o 2JM , , ^, 1 , 

where r' = a/x^ + y'^ + z'^. This /3 is involved in all the dispersion relations. 
As the function /3 changes, the value of the dispersion relation (i.e., the 
relationship between the wave number and angular frequency) changes. In 
this work, we consider only the propagation which changes due to change in 
z and hence the change in /3. 

According to these assumptions, the MHD flow is intended to model a 
wind flowing out of the ergosphere of a Kerr black hole and also an accretion 
flow down to the horizon, both originated at the pair production region, i.e., 
2; = 0. At the pair production region, /3 = — 1. The renormalized shift vector 
is C = P + Vf where Vp = 1, is a flow constant. There is no flow at 2; = 0, 
thus C = I3 + Vf = -I + 1 = at z = 0. 

Our assumed plasma's magnetic field having a variable x-component XB 
with variable A depends on the variable z. The plasma is frozen-in to the in- 
field and we have assumed an adiabatic flow. For simplification, we have also 
taken A = 1, i.e., the x-component of the magnetic field becomes constant and 
leads to a constant unperturbed magnetic field. When we apply these values 
to Eg. (13. 3 p and then substitute the resulting value of V in Eq. (l3.4p with the 
simplification A/po = 1 (as done in [34j-i36j), we obtain the following pairs 
of velocity components. 

1 12 1 
Ml = — tanh(z) /, Vi = -tanh(z) /, (3.5) 

1 12 1 

U2 = -- tanh(z) + -/, V2 = - tanh(2;) + -/, (3.6) 



where / = \/3 — 2tanh (z). The corresponding Poynting vector, then, takes 
the form 

S = — E X B = 2tanh(2)(ex-ez). (3.7) 

An 
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These assumptions lead to a rotating black hole magnetospheric fluid which 
admits a constant external magnetic field. The fiuid is moving with a velocity 
V in the magnetosphere. The gravitomagnetic waves and pair of particles 
are produced in the z = region. The energy extraction from the black hole 
is possible if the medium allows the waves and particles to pass through it 
and move out of the magnetosphere. The gravitomagnetic waves transform 
information within plasma. If the medium, in the region under consideration, 
allows the waves to move normally, these waves can get out of the magneto- 
sphere. This can be well understood by investigating properties of the waves 
in this region. 

The plasma fiow in the magnetosphere is perturbed due to gravity and 
rotation of the black hole. This fiow is characterized by fiuid's density p, 
velocity V and magnetic field B. The linear perturbations will be in xz- 
plane and thus are dependent on x, z and t. The perturbed variables take 
the following form 

P = P^ + pP, V = V° + v, B = B° + 5b, (3.8) 

where unperturbed quantities are denoted by the superscript zero. Further, 
the dimensionless notations p, v and b can be expressed as 

p = — = p{t,x,z), 

P 
V = 6Y = v^(t,x,z)e^ + v^(t,x,z)e:^, 

h = — = h^{t,x,z)e^ + b^{t,x,z)e^. (3.9) 

The perturbations are assumed to be harmonic, i.e., of the form S = 

^-i(ut-kxx-k^z) i^Qy^Yig sinusoidal dependence on t, x and z). 

p(t,x,z) = Ci6, v^(t,x,z) = C2S, Vz(t,x,z) = C36, 

b^(t, X, z) = C46, b;,(t, X, z) = C5S, (3.10) 

where Ci, C2, C3, C4, C5 are arbitrary constants, k^ and kz are the x and 
z-components of the wave vector k = {k^, 0, kz) and u is angular frequency 
of the wave. 



4 Fourier Analyzed Perturbed GRMHD Equa- 
tions 

The perfect GRMHD equations for the Kerr planar analogue (Eq.f l2.2p ) is 
given by Eqs.(2.4)-(2.8) of ^SSj. Substituting the value of rest-mass density 
from Eq.(3.1), we obtain their specific form for the cold plasma. The insertion 
of perturbed values from Eq. fl3.8p and then from Eq. fl3.9p give the component 
form of these equations. This is the same procedure as given in [19] and [3^ . 
Since the perturbed variables have harmonic dependence on x, z and t, 
hence Eq. fl3.10p is used and the following Fourier analyzed form is obtained. 

LkzC2 — {ik-^X + A')c3 — Ci^tkzU — lu + u') + c^irriLkz + m!) = 0, (4.1) 
ikxC2 — ikx\c?j + Lc^{mkx + kzU — tj} = 0, (4.2) 

kxC4: = -kzC5, (4.3) 

Cil{—uj + mkx + ukz\ + C2 \—lw'~^'^V + ik^j'^uV + ikxq+ 
7^m{(1 + 27V2)\/' + 2^'^uVu'}] + C3 [~l{uj + pk)^'^u 



(4.4) 



+Lk,r + LkxYuV + 27Vl/\/' - (1 - 27^M^)r- 

u 

Cip'j'^u{qV' + 'j'^uVu'} + Ci \_—iw {p'y'^q + s^ + cy [p'j'^q — s} 

+P7^m{(1 + 4-f^V^)uu + AqVV'}] + C3 [-iw {p-f\V - As} 

+Ly {p-f\V + As} + P7'[{(1 + 272m2)(1 + 27V) 

-7V2}V' + 27^(1 + 2j\^)uVu'] + uX's] + sci{-tk,{l - u^) 

+uu'} + SC5 {-A' - um' + Lkx{l - V^) - 2iuVk,} = 0, (4.5) 

ci'j'^ p[u{ru' + 'j'^VuV'} — V(3'] + C2 [—iw {p'j'^uV — As} 

+Ly {p^%V + As} + p72{7VV^'(l + A^^V^) - /3'{1 + 27^^^) 

+2V-f^uu'{l + 27^^^)}] + C3 [-L {p7^(l + 7^^^) + A^s} w 

+Ly {p7V - A^s} + P7^[m'(1 + 7^m^)(1 + 47^m^) 

+2uVj^{m' + 2-f'^u^V'}] - sWu] + SC4 {iA;^A(l - u^) 

+A' - Xuu'} + sc5{2XuVLk, + Xum - tk^Xil - V^)] = 0, (4.6) 



where m = V — (3, q = 1 + 7^^^^, r = 1 + 7^^^, s = ^, w = ijJ + (5kx and y = 
k^V + k^u. Equation (14. 3 p gives the relation between k^ and k^, i.e., k^ = 
— f-kz- We assume C4 = C5 yielding kx = —kz which reduces our wave vector 
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to k = (kx, 0, — A;^) with wave number A; = |k| = A/2AJf. This assumption will 
make the set of equations easier to solve. The relaxation in this assumption 
may change the regions of normal and anomalous dispersion. 

5 Numerical Solutions 

When we solve the determinant of the coefficients of Eqs.(4.1), (14. 2p . (14. 4p - 
(14. 6 p of constants (ci, C2, C3, C4, C5) for the x-component of the wave vector, 
it gives a complex dispersion relation of the form: 

Ai{u, z)k/ + A2{uj, z)kj' + A^{u, z)kj + ^4(0;, z)k^ + A^iu, z) 
+i{Bi{uj, z)kj' + 52(w, z)k^^ + B-i{uj, z)kj^ + ^4(0;, z)k^^ 
+B^{uj,z)k^ + BQ{u,z)} = Q. (5.1) 

This relation is quintic in kx and cannot give exact solutions. Thus, we 
solve it numerically by using Mathematica, for the complex values of kx- The 
assumption kx = —kz implies that kz is also a complex number. Thus the 
quantities given by Eq. (l3.10p take the form 

where kx = kjix + tkjx and kz = kjiz + ikjz- The values k^x and k^z represent 
the X and z-components of the propagation factor from which we can obtain 
the X and z-components of the phase and group velocities whereas kjx and 
kiz represent the x and z-components of the attenuation factor. 

We have considered the region —5 < 2; < 5 to analyze the wave properties 
of the cold plasma. The cold flow assumption indicates that the fluid has no 
viscosity and no heat conduction. This is indeed an ideal case but provides 
basis for the case of hot plasma, thus has its own importance. Perturbations 
themselves modify this assumption because we include the perturbations in 
the GRMHD equations for cold plasma. Since the flow variables m, V, A 
have more variations near the pair production region, thus we have excluded 
the region — 1 < 2; < 1 and solved the relation for rest of the region where 
these variables have small variation. Thus the region under consideration 
—5 < z < 5 is divided into two regions —5 < z < —1 and 1 < z < 5 for 
which roots are numerically interpolated. The region —5 < z < —1 indicates 
the neighborhood of the pair production region towards the event horizon 
and the region 1 < 2; < 5 shows the neighborhood of the pair production 
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region towards the outer end of the magnetosphere. We take the step-length 
0.2 for z and u and find the value of k^ at each point (zj, Wj) (i = 1, ..., 21, j = 
1, ...51) of the mesh. Then by separating each root, we estimate the surface 
by numerical interpolation. The a;-components of the propagation vector, 
attenuation vector, phase and group velocities can be evaluated from these 
interpolation functions. 

We have obtained two velocities (given by Eqs. fl3.5p and (13. 6p ) for our 
assumed plasma (i.e., rotating cold plasma with constant rest-mass density). 
For each velocity, Eq. (l5.ip leads to five values of the x-component of the wave 
numbers. We display x-components of the propagation vector, attenuation 
vector, phase and group velocities and analyze the dispersion with respect to 
these quantities. 

For the region towards the black hole horizon with the fluid flow velocity 
given by Eq. (l3.5p . we obtain one root for which the wave number is infinite 
at several values of z. Such a root is also found for the region towards 
the outer end of the magnetosphere [z = oo) for the velocity components 
(Eq. (l3.6p ). These roots indicate that the waves are evanescent there which 
decreases with the passage of time. Thus the energy cannot pass through 
the region via waves and hence these cases are not interesting. For the 
velocity components (13. 5p . the dispersion relations obtained in the region 
towards the event horizon in the neighborhood of the pair production region 
are given by Figures 1-4. Figures 6-9 indicate the dispersion relations for the 
region towards the outer end of the magnetosphere. The dispersion relations 
related to the velocity components (13. 6 p for the region towards the event 
horizon are given by Figures 10-14. For the region towards the outer end of 
the magnetosphere, these relations are shown by Figures 15-18. 

In Figures 1-4 and 10-14, the positive propagation factor shows that the 
waves are moving towards the pair production region whereas the negative 
propagation factor indicates that the waves are moving towards the event 
horizon. We see from Figures 6-9 and 15-18 that when the propagation factor 
is positive, the waves are moving towards the outer end of the magnetosphere 
whereas the negative propagation factor indicates that the waves are moving 
towards the pair production region. Table 1 indicates the regions of magne- 
tosphere with the directed waves. The rest of the region contains random 
points of positive and negative propagation factor and thus the direction of 
waves changes at each point. 

If the attenuation factor is increasing, it leads to the damping of waves 
whereas if this factor decreases, it indicates wave growth. Table 2 shows 
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Figure 

No. 


Region of waves 

moving towards 

the event horizon 


Region of waves 

moving towards 

pair production region 


Region of waves 

moving towards 

the outer end 


1 


— 


-2.62 < z < -1* 


— 


2 


-2 < ^ < -1 


— 


— 


3 


— 


-2< ^< -1 


— 


4 


— 


-2.1 <z<-l 


— 


5 


— 


Entire Region 


— 


6 


— 


— 


Entire Region* 


7 


— 


— 


Entire Region* 


8 


— 


— 


Entire Region* 


9 


— 


— 


Entire Region* 


10 


Entire Region 


— 


— 


11 


— 


Entire Region 


— 


12 


— 


Entire Region* 


— 


13 


— 


Entire Region* 


— 


14 


— 


Entire Region 


— 


15 


— 


1 < z < 2.7 


— 


16 


— 


— 


l<z< 2.25* 


17 


— 


l<z<2* 


— 


18 


— 


— 


Entire Region 



Table 1: Table indicating the direction of waves in respective regions. 



the damping and growth of waves. The up-arrow in the table shows the 
increase whereas the down-arrow shows the decrease in respective quantity. 
The attenuation factor takes random values otherwise. 

In a region where the phase velocity is greater than the group velocity, 
the waves are dispersed normally [5S]-[Sn] whereas the group velocity greater 
than the phase velocity shows that the waves disperse anomalously. Table 3 
shows the regions of normal and anomalous dispersion of waves in respective 
figures. Rest of the region contains random points of normal and anomalous 
dispersion. 

We note that in all the tables star (*) means validity of the attribute to 
the whole region except for waves with negligible angular frequency. 
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Figure 

No. 


Regions of 
Wave Growth 


Regions of 
Wave Damping 


1 


-1.6 <z< -1,* 

0.398 <u< 7.25 

as w t ci'iid 2; 4, 


— 


2& 3 


Random 


Random 


4 


— 


-1.5 < 2; < -1 as cj t and z | 


5 


< w < 0.23 as w t 
1 < 2< 1.4, 0.23 < w < 10 

as 2; t 


1.4 <z < 5 aszt 
l<z< 1.4, < w < 0.23 

as 2; t 


6 




1 < ^ < 2, 0.35 < w < 10 

as 2; t and uj \. 

2 < z < 5, 0.35 < w < 10 

as z t 


^7 

i 


— 


As ^t* 


8 


— 


As a; t and z \. 


9 


1.425 < ^ < 5, 0.2 < cj < 10 

as 2; t 
1 <z <5, < w < 0.2 

as 2; t 


l<z< 1.425, 0.2 < cj < 10 
as 2; t 


10 


-1.35 < z < -1 as ^ t 
1 < 2 < 2, 0.195 <tj < 10 

as 2; t 


-5 < 2 < -1.35 as z\ t 


11 


-1.35 <z<-l as zi 


-5 < z < -1.35 as ^i 


12 


As \z t 


As w t 


13 


z \. and 00 t 


— 


14 




-h<z< -1,0.175 < w < 10 

as 2; t and uj t 


15 


1.6 < 2< 1.8, <cj < 0.195 
as z t and 00 t 


1 < 2< 2, 0.195 < w < 10 

as 2; t 
l<z< 1.6, < w < 0.195 

as 2; t and u t 
1.8 < ^< 2, < w < 0.195 

as z t and u t 


16 


1.5 < 2 < 1.625 


1.0<2< 1.5 


17 & 18 


Random 


Random 



Table 2: Table indicating the regions of growth and damping of waves. 
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Figure 

No. 


Regions of 
Normal Dispersion 


Regions of 
Anomalous Dispersion 


1 


-1.6 < z < -1, 0.1 <a; < 7.25 


— 


2 


— 


-2.15 < 2 < -1, <a; < 0.05 


3 


Random 


Random 


4 


— 


-2 < z < -1 


5 


l<z< 1.25, <cj < 2 
1.25 <^<3, 0<cj<l 
3<2<5, 0<a;< 0.75 


— 


6 


— 


2< ^ < 5, 0.8 < w < 1, 

1 < 2 <4, 0.5 < tj < 0.8, 

1 < 2 < 5, <a; < 0.5 


7 


Entire region* 


— 


8 


— 


0.2225 < w < 5 


9 


uj < 0.4025 


— 


10 


-b<z< -4.2, <a; < 0.1 


— 


11 


-3 < ^ < -2.5, < cj < 0.4 


-5 < ^ < -3, < w < 5 
-2.5 < 2 < -1, <cj < 10 


12 


Random 


Random 


13 


-5 < z < -2, < w < 0.2 


-5 < 2< -1, 0.2 < w < 10 
-2<z< -1, 0.16 <a; < 0.2 


14 


-2<z< -1, 0.1 <cj < 0.39 
-2<z< -1, 0.41 < tj < 0.5 
-2.2 < 2 < -1, <cj < 0.08 


— 


15 


— 


1.8 < ^ < 2, 0.002 <uj< 0.003 


16 


1 < 2 < 1.4, 2.5 <w < 10 
l<z< 1.4, 1 < w < 1.6 
l<z< 1.25, 1.6<a; < 2 


— 


17 


1.8 < 2 < 2, < w < 0.00015 


l<z< 1.5, 0<u< 1.25 


18 


1 < 2 < 1.45, < cj < 3.5 


— 



Table 3: Table indicating the regions admitting normal and anomalous dis- 
persion of waves. 
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6 Conclusion 

This work is devoted to investigate the wave properties in the neighborhood 
of the pair production region. We have considered rotating cold plasma filled 
magnetosphere. The perturbations are assumed to be simple harmonic waves 
caused by the gravity infiuence as well as the rotation of the black hole. The 
determinant of the Fourier analyzed GRMHD equations is solved to obtain 
the x-component of the wave vector which is found to be a complex number. 
This has been done numerically for the two values of the velocity of the 
magneto-rotational fiuid obtained in Section 3. The relation between the x 
and 2;-components of the wave vector is k^ = —k^. This indicates that if we 
can have information about k^ and its corresponding quantities, we can infer 
the results for the waves moving along 2;-axis. 

The summary of these results can be expressed as follows: 

The propagation and attenuation factors take random values far from 
the pair production region in Figures 1-4, 15-18. The propagation factor 
increases with increasing angular frequency in Figures 6-9, 11-14. This gives 
the increment in propagation of waves along x-axis with the increment in 
angular frequency. It is observed in Figures 5-14 that the phase and group 
velocities take their extreme values (either maximum or minimum) near the 
pair production region. Thus the pair production region allows the waves to 
take extreme values near it. 

In most of the figures, the region in the extreme neighborhood of the 
pair production region admits the waves to pass through (normal dispersion) 
which indicates that the waves can pass through this region. Anomalous dis- 
persion in the extreme neighborhood of the pair production region in Figures 
2, 4, 6, 8, 11, 13 and 17. In the region of anomalous dispersion, the phase 
velocity is less than the group velocity, so the waves cannot carry energy 
alongwith. Most of the figures show random dispersion of waves in the far 
regions of the pair production region. 

Figure 7 shows normal dispersion of waves throughout the region which 
indicates that the waves can move towards the outer end of the magneto- 
sphere. This verifies the result given by Zhang [33] and Sheikh |37] that cold 
plasma allows outfiow of energy fiux from the pair production region. 

It would be worth exploring to check these properties by taking a hot 
plasma in a rotating magnetosphere. 
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Figure 1: A large region admits normal dispersion of waves moving towards 
the pair production region. 
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Figure 2: A small region in the neighborhood of the pair production region 
indicates anomalous dispersion at small angular frequencies. A small region 
shows the waves are moving towards the event horizon. 
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Figure 3: Random dispersion of waves is found. A small region shows that 
the waves are moving towards the pair production region. 
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Figure 4: The region near to the pair production region shows anomalous 
dispersion for the waves moving towards the pair production region. 
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Figure 5: Wave propagation decreases as the waves move towards the pair 
production region. Small regions admitting normal dispersion are found. 
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Figure 6: The waves damp as they move away from the pair production re- 
gion except for the waves with very small angular frequencies. Regions of 
anomalous dispersion are found for small angular frequencies in the neigh- 
borhood of the pair production region. 
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Figure 7: The waves damp while moving away from the pair production re- 
gion (except for the waves with small angular frequencies). Normal dispersion 
is found except for the negligible angular frequency waves. 



19 



Propagation Factor 



Attenuation Factor 




Phase Velocity 



Group Velocity 



lit 0.999 
0.99. 
0.998 





Figure 8: Waves grow as they move away from the pair production region. 
In most of the region, dispersion is anomalous except for the waves with very 
small angular frequencies lying between 0.2225 and 5. 
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Figure 9: The propagation of waves increases with increasing angular fre- 
quency. Waves with small angular frequencies admit normal dispersion. 
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Figure 10: Waves grow and then damp when the value of z increases. This 
propagation decreases as their angular frequency increases. Dispersion is 
found to be normal in a small region. 
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Figure 11: The propagation of waves decreases when the waves move away 
from the event horizon. Waves damp first and then grow while moving to- 
wards the pair production region. Small regions show normal and anomalous 
dispersion of waves. 
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Figure 12: Waves damp as they move towards the pair production region. 
Random points of normal and anomalous dispersion are found. 
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Figure 13: The increase in the angular frequency increases the wave prop- 
agation. Waves grow as they move towards the pair production region. A 
region of normal dispersion is found for small frequency waves. 
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Figure 14: The wave propagation decreases and the waves damp as they 
move towards the pair production region except for the low angular frequency 
waves. Small region admitting normal dispersion of waves are identified. 
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Figure 15: The medium admits a small region for the waves moving towards 
the pair production region. The wave propagation decreases and then in- 
creases as the waves move away from the event horizon. Small regions of 
normal and anomalous dispersion are found. 
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Figure 16: The propagation vector, attenuation factor, phase and group 
velocity vectors admit random values away from the pair production region. 
Small regions near the pair production region show normal dispersion. 
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Figure 17: A small region admitting anomalous dispersion of waves moving 
towards the outer end of the magnetosphere. 
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Figure 18: A small region in the neighborhood of the pair production region 
shows normal dispersion of the low angular frequency waves moving towards 
the outer end of the magnetosphere. 
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